Slave-rotor mean field theories of strongly correlated systems 
and the Mott transition in finite dimensions 
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The multiorbital Hubbard model is expressed in terms of quantum phase variables ( "slave rotors" ) 
conjugate to the local charge, and of auxiliary fermions, providing an economical representation of 
the Hilbert space of strongly correlated systems. When the phase variables are treated in a lo- 
cal mean-field manner, similar results to the dynamical mean-field theory are obtained, namely a 
Brinkman-Rice transition at commensurate fillings together with a "preformed" Mott gap in the 
single-particle density of states. The slave- rotor formalism allows to go beyond the local description 
and take into account spatial correlations, following an analogy to the superfluid-insulator transition 
of bosonic systems. We find that the divergence of the effective mass at the metal- insulator transi- 
tion is suppressed by short range magnetic correlations in finite- dimensional systems. Furthermore, 
the strict separation of energy scales between the Fermi- liquid coherence scale and the Mott gap, 
found in the local picture, holds only approximately in finite dimensions, due to the existence of 
low-energy collective modes related to zero-sound. 
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I. INTRODUCTION 

Strongly correlated fermion systems constitute a chal- 
lenge, both from a fundamental point of view (with phe- 
nomena such as the Mott transition [lj and high- tem- 
perature superconductivity), and on a more quantitative 
level with the need of reliable tools to handle intermediate 
and strong coupling regimes (even for simplified models 
such as the Hubbard model). In recent years, the dynam- 
ical mean-field theory (DMFT) has allowed for significant 
progress in this respect |2|. In particular, this approach 
has led to a detailed theory of the Mott transition, and to 
a quantitative description of the physics of strongly cor- 
related metals. Despite these successes, the limitations 
of this approach have been emphasized on many occa- 
sions. The main one has to do with the effect of spatial 
correlations (e.g magnetic short-range correlations), and 
more precisely with the effect of these correlations on the 
properties of quasiparticles. For example, the tendency 
to form singlet bonds due to superexchange is widely be- 
lieved to be a key physical effect in weakly doped Mott 
insulators. Also, at the technical level, the application 
of DMFT to materials with a large orbital degeneracy 
(e.g in combination with ab-initio methods 0,0, 01), as 
well as cluster extensions of DMFT 0, El are computa- 
tionally challenging because they involve the solution of 
a multi-orbital quantum impurity model. 

For these reasons, there is still a strong need for 
approximate, simpler treatments of strongly correlated 
fermion models. Those treatments should incorporate 
some of the DMFT successes (e.g regarding the descrip- 
tion of the metal- insulator transition), but they should 
also pave the road for describing physical effects beyond 
DMFT at least at a qualitative level. 



The purpose of this paper is to present a simple mean 
field description of correlated systems which fullfills some 
of these goals. Our main idea is to focus on the degrees of 
freedom associated to the relevant physical variable as- 
sociated to the Mott transition, namely a slave quantum 
rotor field, dual to the local electronic charge. This slave 
rotor representation was introduced previously by us for 
the description of quantum impurity models and meso- 
scopic devices 0, 0) an d is applied here in the context 
of lattice models. This allows for a simple reformula- 
tion of the orbitally degenerate Hubbard model, which, 
when the interaction has full orbital symmetry, is quite 
superior to previously developed slave- boson represen- 
tations HEaim. 

When the simplest possible (single-site) mean-field ap- 
proximation is used in conjunction with this slave-rotor 
representation, a description of the Mott transition very 
similar to that of DMFT is found. The metallic phase 
disappears through a Brinkman-Rice transition, at which 
the quasiparticle wieight vanishes and the effective mass 
diverges. The slave rotor approach does preserve Hub- 
bard bands in the insulator, and a "preformed" Mott 
spectral gap opening up discontinuously at the transition 
is found, as in DMFT. 

The most interesting aspect of our approach lies how- 
ever in the possibility of going beyond this purely lo- 
cal mean-field description. By decoupling the spinons 
and slave rotor degrees of freedom, the Hubbard 
model is mapped onto a free spinon hamiltonian self- 
consistently coupled to a quantum XY lattice model. The 
(dis)ordering transition of the latter corresponds to the 
Mott transition, in analogy with the superfluid-insulator 
transition of the bosonic Hubbard model. Because spatial 
correlations are now included, we find important modifi- 
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cations to the DMFT picture. In particular, the effective 
mass remains finite at the transition, due to the quench- 
ing of the macroscopic entropy by magnetic correlations 
in the Mott phase. Importantly, low-energy charge col- 
lective modes are shown to affect the opening of the Mott 
gap, which now develops in a continuous manner, so that 
the separation of energy scales found in DMFT only holds 
in an approximate manner. These simple results can 
be considered as deviations from the DMFT predictions 
that could possibly be observed in photoemission exper- 
iments 0]. However, restauration of the local gauge 
symmetry should occur due to fluctuations beyond the 
mean-field approximation, possibly modifying the latter 
result on a qualitative way. 

The paper is organised as follows. In section [H] we 
introduce the exact slave rotor description of a simple 
atomic level with orbital degeneracy, and show that an 
approximate treatment of the local constraint is sufficient 
to describe correctly the full Coulomb staircase, as well 
as one- particle spectra. Then, in section Hill we develop 
the simplest (local) mean field treatment of both the An- 
derson and Hubbard models, and in the latter case, study 
the multiorbital Mott transition. Finally, spatial fluctu- 
ations beyond DMFT are included in section IIVI with 
an emphasis on the behavior of the effective mass and 
the excitation spectrum. The conclusion presents sev- 
eral possible applications and extensions of our formal- 
ism, and also discusses some of the open issues raised by 
our results. 



II. ROTOR REPRESENTATION OF 
INTERACTING FERMIONS 

A. Slave-rotor representation 

In Ref. (see also 0), we introduced a represen- 
tation of the Hilbert space of N fermions d\, in terms 
of a collective phase degree of freedom 0, conjugate to 
the total charge, and of N auxiliary fermions /J. The 
spin/orbital index runs over N values a — 1 ... N (e.g 
a =t, | for N — 2). In the following, we consider only 
interactions which have the full SU (N) spin/orbital sym- 
metry. Let us consider the Hamiltonian corresponding to 
a single "atomic level" , in the presence of a local Hubbard 
repulsion: 



(i) 



The crucial point is that the spectrum of the atomic 
Hamiltonian Q depends only on the total fermionic 
charge Q = 0, • • • , N and has a simple quadratic depen- 
dence on Q: 



Eq =e a Q+j 



(2) 



There are 2 N states, but only N + 1 different energy lev- 
els, with degeneracies (g). In conventional slave boson 
methods 0, 0], a bosonic field is introduced for each 
atomic state \a\ ■ ■ ■ <jq) (along with spin-carrying aux- 
iliary fermions /J). Hence, these methods are not de- 
scribing the atomic spectrum in a very economical man- 
ner, and lead to very tedious calculations when orbital 
degeneracy becomes large, even at the mean-field level 

Ea m m 

However, we stress that the extensive 
number of degrees of freedom necessary in those other 
approaches can become useful when the SU (N) symme- 
try is broken, either by magnetic order or crystal fields. 

The spectrum of Q can actually be reproduced by 
introducing, besides the set of auxiliary fermions /J, a 
single additional variable, namely the angular momen- 
tum L = —id /dO associated with a quantum 0(2) rotor 
9, an angular variable in [0, 2n]. Indeed, the energy levels 
@ can be obtained using the following Hamiltonian 



(3) 



A constraint must be imposed, which insures that the 
total number of fermions is equal to the 0{2) angular 
momentum (up to a shift, in our conventions): 



E 



fafcr 



(4) 



This restricts the allowed values of the angular momen- 
tum to be £ = Q - N/2 = -N/2, -N/2 + 1, • • • , N/2 - 
1, N/2, while in the absence of any constraint £ can be an 
arbitrary (positive or negative) integer. The spectrum of 
© is e Q + U£ 2 /2, with £ = Q- N/2 thanks to Q, so 
that it coincides with 

It is easily checked that the full Hilbert space is cor- 
rectly described as: 



\a 1 ...a Q ) d = \a 1 ...a Q ) f \£ = Q-N/2) t 



(5) 



in which \a± . . . <jq) d t denotes the antisymmetric fermion 
state built out of d- and /-fermions, respectively, and 
\£}g denotes the quantum rotor eigenstate with angu- 
lar momentum I, i.e. (9\£) g — e ae . For N = 2, 
this corresponds to: = |T) / |0) e , = \i) f \0) g , 

\U) d - |U)/|+l> fl and |0) d = 10^1-lV The creation 
of a physical electron with spin a is associated to the 
action of fl on such a state as well as raising the to- 
tal charge (angular momentum) by one unit. Since the 
raising operator is e lB , this leads to the representation: 



4 = fU 



-iO 



(6) 



The key advantage of the quantum rotor representation 
is that the original quartic interaction between fermions 
has been replaced in J3Jl by a simple kinetic term for the 
phase field, (U/2)L 2 . 

We point out here that a similar phase representation 
was developed before in the context of Coulomb block- 
ade in mesoscopic systems, see e.g 0,0,0]- However, 
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the present work and our previous paper present the 
first applications of the rotor technique to the context 
of strongly correlated lattice models. In particular, the 
question of quasiparticle coherence which is crucial to the 
description of a Fermi liquid cannot be investigated seri- 
ously with a phase-only description , as shown in Q . 
In this perspective, the slave rotor should be seen as a 
natural extension (and simplification) of the usual slave 
boson techniques |20l l2l| in the context of a finite but 
orbitally symmetric Coulomb repulsion. In principle, it 
can also be applied to systems with long-range interac- 
tions mm. 



B. Treating the constraint on average: atomic limit 

In the following, we will study different kinds of mean- 
field approximations based on this slave-rotor represen- 
tation. A common trait of these mean-field approxima- 
tions is that the number constraint (0J will be treated 
on average. This is equivalent to treating the constraint 
in a "grand-canonical" ensemble, which would of course 
be exact in the limit of a large spin/orbital degeneracy 
N — > oo. In this section, we investigate the accuracy of 
this approximation for the atomic Hamiltonian IjTJl. for 
finite values of N. 



1. Coulomb staircase: occupancy vs. eo 

Let us first consider the dependence of the average oc- 
cupancy (Q) — (X)cr d\d a ) on the position of the atomic 
level eo, which reads: 



(Q)at 



1 N 

— y 



2=o 



Qe 



-0E Q 



(7) 



with Z a t = 2^2q (cj) e ^ Eq ■ In the limit of zero tempera- 
ture, the dependence of (Q) on eo is the "Coulomb stair- 
case" in Fig.^ When treating the constraint on average, 



-3U/2 -U/2 



(Q) 



at 



U/2 3U/2 



eo 



FIG. 1: Coulomb staircase in the atomic limit for the case of 
two orbit als, N — 4. 

a Lagrange multiplier h is introduced which is conjugate 
to J3J, and one optimizes over h instead of fully inte- 
grating over it. This amounts to consider the following 



effective Hamiltonians: 



HS 



2 



,.u-f 



(8) 



(9) 



The Lagrange multiplier h is determined by the average 
constraint equation: 



(L) h = N 



n F (e ~h)-~ 



(10) 



in which (L)h is the average of L in the Hamiltonian 10: 

^9 , — ' 



with: E t = U£ 2 /2 + hl and Zg = J2e e ^ PEe - Solving (HJ 
for ft as a function of eo and temperature T = 1/(3 yields 
the dependence of the total charge within this approxi- 
mation: 



(Q) =Nn F [e -h(e ,T)} 



(12) 



We need to compare this approximation to the exact re- 
sult Q in the atomic limit. A graphical representation 
(Fig. |2J is useful in order to understand the solution of 
(|TUj) . At T = 0, one finds h = e , as long as < Q < N. 



inp(n) — h) — 2 



J0_ 



Up, ' 3*7/2 



-3*7/2 -U/2 
-1 

-3 



FIG. 2: Graphical solution of the average constraint equa- 
tion HUH . The intersect (cross) moves exactly along the 
Coulomb staircase shown in figure 

The exact dependence of the average charge Q upon eo 
at T — is correctly reproduced by our approximation, 
corresponding to the "Coulomb staircase" : 



N 
~2 



= — + I for -U 



21+1 TT 21 - 1 f^S 

— — < e < -U — — (13) 



Note that h — eo vanishes linearly with temperature ac- 
cording to: h = eo — Tln(N/Q — 1) + • • • : this is why 
the full Coulomb staircase can be reproduced with a sin- 
gle Fermi factor in At finite temperature, our ap- 
proximation does not coincide with the exact result for 
Qat (eo ,T), but deviations are only sizeable for tempera- 
tures comparable to U , which is not a severe limitation 
in practice. 
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2. Spectral functions 

We now study the consequences of the approximate 
treatment of the constraint for the Green's function and 
spectral function. Following Q8l9p . the quantum rotor and 
auxiliary fermion degrees of freedom are described by two 
independent Hamiltonians, so that the Green's function 
of the physical electron Grf(r) = — (Td (T (r)d^(0)) factor- 
izes into: 

G d {T) = Gf{T)G e {r) (14) 

with G 6 {t) = (expi [6(0) - 0(t)}). Equivalently, the 
physical electron spectral function is given by: 

Pd{u) = - J duj' pf(w')pe(uj-uj') [tif(—u)') + n s (w - u')] 

Let us consider T = 0, and eo in the range corresponding 
to the plateau of charge Q in the Coulomb staircase. The 
ground-state energy is Eq = U(Q — N/2) 2 + e^Q and its 
degeneracy is (q). The two excited states obtained by 
adding or removing a particle correspond to transition 
energies: A± = E Q±1 - E Q = ±e ± U{Q - N/2 ± 1/2). 
When acting with d\ on the ground-state, only those 
ground-state components which do not already contain 
a contribute, and there are i^q 1 ) = (1 — n){q) sucn 
components. Similarly, when acting with d a , only the 
components in which a is occupied contribute, and there 
are (q~i) = n (q) °^ them. From these considerations, 
we see that the exact spectral function reads, at T = 0: 

pf(u) = (l - |) % - A+) + | 6(w + A_) (16) 

These two atomic transitions are the precursors of the 
Hubbard bands in the solid. Note that they have unequal 
weights, except at half-filling Q — N/2. At finite tem- 
perature, additional peaks appear (except for N = 2), 
corresponding to transition between two excited states 
(with exponentially small weight for T -C £7) . 

Remarkably, the expressions I|14ll5fl in which the quan- 
tum rotor and auxiliary fermions are treated as decou- 
pled, do reproduce this exact result at T = 0. The 
easiest way to see this is to notice that, at T = 0, 
Gf(r) = -(1 - Q/N)8(t) + (Q/N)8(-t), since n F {e - 
h) = Q/N. The rotor Green's function G 9 is e~ A+T 
for r > and e A - T for r < 0. Substituting into {2}, 
this corresponds to the exact expression i|16[l . Alterna- 
tively, one can use the expressions of the T = spec- 
tral functions into (|15fl : pf(w) — S(lo — eg + h) and 
po{w) — —8(oj — A_|_) + 8(u> — A_), keeping in mind that 
^f(co — h) = Q/N while jib(w — eo + h) = —9(—lu) 
as T — » 0. Again, deviations between the approximate 
treatment and the exact results are found at finite tem- 
perature, but remain small for T <S U. Let us emphasize 
that, because the rotor Green's function Gg is a continu- 
ous function at r = 0, with Gq(t = 0) = 1, the factorized 
approximation l|14|) insures that the physical (d-electron) 



spectral function is correctly normalized with total spec- 
tral weight equal to unity. 

To summarize, we have found that treating the con- 
straint on average reproduces accurately the atomic limit 
at T = 0, both regarding the Coulomb staircase depen- 
dence of Q vs. eo, and regarding the spectral function. 
This is a key point for the methods introduced in this ar- 
ticle, which allows them to describe reasonably the high 
energy features of strongly correlated systems. 



3. Functional integral formulation 

We briefly introduce here a functional integral formal- 
ism for the fl and 8 degrees of freedom, and derive the 
action associated with . This is simply done by switch- 
ing from phase and angular-momentum operators (8, L) 
to fields (8, d T 8) depending on imaginary time r € [0, 0\, 
with 9(/3) = 9(0) + 2irn. The action is constructed from 
S = j^dr[—iL d T 8 + H + pd T f] , and an integration over 
L is performed. It is also necessary to introduce a La- 
grange multiplier h in order to implement the constraint 
L = J2 a f\fo — N/2. We note that, because of the charge 
conservation on the local impurity, h can be chosen to be 
independent of time, with h £ [0,2tt/{3]. This leads to 
the following expression of the action: 

5 at = jy^m +e + ih)U + {dr9 2 + h)2 - ijh 

The constraint is implemented exactly provided h is inte- 
grated over. The above approximation amounts to eval- 
uate the integral by a saddle-point approximation over 
h, and the saddle-point is found to be on the real axis, 
with h = ih. 

Finally, let us mention that, in a previous publi- 
cation [3, we have explained in detail the connec- 
tion between the rotor construction and the Hubbard- 
Stratonovich decoupling of the interaction in the charge 
channel. 



III. THE SIMPLEST MEAN-FIELD 
APPROXIMATION 

In this section, we introduce a very simple mean-field 
approximation based on the slave rotors variables. This 
approximation is similar in spirit to the condensation of 
slave bosons in conventional slave boson mean-field the- 
ories. We illustrate this approximation on two examples: 
the Anderson impurity model and the Hubbard model. 
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A. Anderson impurity model 

The Anderson impurity model describes a local orbital 
hybridized to a conduction electron bath: 

H = H at + J2 aicte + V 5>LA + h-c) (17) 



ka 



ka 



This Hamiltonian can be rewritten in terms of the slave 
rotor and auxiliary fermion variables: 

a ka 

+ Vj2(clUe-* e + h.c.) (18) 

ka 

submitted to the constraint (@J. The simplest possible 
approximation is to decouple the rotor and fermion vari- 
ables, leading to two effective Hamiltonians: 

Hf = ( eo - h )J2fafa+J2^ C l C «° 

(7 ka 

+ ^eff 5>L/a +h.c) (19) 



ka 



U 



Hf = ^-L 2 + hL + K cos( 



(20) 



The parameters V e g, K and h in these expressions are 
determined by the coupled self-consistent equations: 



V cS = V(cost 



ka 



(L) e = N 



n F (e -h)- 



1 



(21) 
(22) 

(23) 



in which the averages are calculated with the effective 
Hamiltonians above. 

Let us first examine the particle-hole symmetric case 
(o = (Q = N/2) in which the solution of (j2HJ is ft, = 0. 
The rotor sector is described by the effective Hamiltonian 
(|20|1 corresponding to the Schrodinger equation: 



2 d9 2 



K cos 9 



*(<?) = EV(9) 



(24) 



For K = 0, the ground-state wave function is the state 
|Z = 0), uniform on [0, 2tt], corresponding to maximal 
phase fluctuations and thus to the absence of charge fluc- 
tuations. This is associated with the atomic limit, as 
explained above. As soon as the hybridization V is non- 
zero, we shall see that K ^ 0. The wave function is 
then maximum (K < 0) for 9 = 0, 2ir, and (cos 9) ac- 
quires a non-zero expectation value. This corresponds to 
a non-zero effective hybridization V c g = V(cos9), so that 
the auxiliary fermion effective Hamiltonian is that of a 



resonant level model. This captures the physics of the 
Kondo effect, and the corresponding Kondo resonance at 
the Fermi level. 

Even though V is a singular perturbation on the atomic 
limit, its effect can be easily understood analytically in 
the present framework by treating the potential energy 
Kcos9 perturbatively. To first order in K, the ground- 
state wave function reads: 



|^ 1) > = |0) 



with (9\£} = e lU and E t 
(co S 9) g = 



-2K 



(£\Kcos9\0) 

En — Et 



10 



(25) 



E 

1^0 



U£ 2 /2. This yields: 

ob el*? 5 ) 

I cos 6> |£) I 2 



Et 



2K 



(26) 



Hence, using (|2ip. one obtains K = —UV c s/2V, which 
yields the following self-consistent equation for the effec- 
tive hybridization V e g when substituted into (|22[l: 



2V 2 
IF 



E 

ka 



{c\ja + flc k<J ) f (27) 



The right-hand side of this equation is easily evaluated for 
the resonant level model l|19|l . For simplicity, we consider 
a flat conduction band £k G [—A, A], and focus on the 
universal regime: A > U > T, with T = -kV 2 /2K. To 
dominant order in 1/A, l|27|l reads: 



1 = —N- 



2V l 



duj- 



UA ./-a w 
2A 2 



(V c %n/(2A)y 



2V 2 
N UA ln 



(28) 



This yields the following expression for V c s and for the 
width of the Kondo resonance when A> U 3> T: 



r, 



eff 



2A 



Aexp 



ttU \ 
4NT ) 



(29) 



This coincides with the exact expression [23]. The local 
orbital spectral function obtained from (|14|) reads: 



r C ff/7r 
^ 2 + r" 



inc 
Pd 



(OJ) 



(30) 



cir 



The first term in this expression is the Kondo resonance, 
and carries a spectral weight Z = T e s/T = (cos 9) e . It 
satisfies the Friedel sum-rule pd(oJ — 0) = I/71T. Away 
from the particle-hole symmetric case (eo 7^ 0), the loca- 
tion of the resonance is set by eo — h, which is the renor- 
malized impurity level familiar from conventional slave- 
boson theories. The rotor approximation does conserve 
total spectral weight, and therefore yields an incoherent 
contribution to the spectral function with a weight 1 — Z. 
This incoherent contribution is correctly centered around 
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the atomic transitions, as explained above. However, the 
width of these Hubbard bands is incorrectly described 
by the simple approximation presented here, in which 
phase fluctuations are underestimated at short times. As 
a result, the Hubbard bands have a bandwidth of or- 
der r e ff in this approximation (instead of the expected, 
and much broader width, of order T). We note how- 
ever that conventional slave-boson approximations with 
a condensed boson neglect altogether the Hubbard bands 
at the saddle-point level, and therefore the present ap- 
proximation, simplified as it may be, is preferable in this 
respect. An improved method for the treatment of phase 
degrees of freedom, leading to a much more accurate de- 
scription of the Hubbard bands, has been discussed in 
previous publications 0, • This method consists in a 
set of coupled integral equations for the Green's func- 
tions of the auxiliary fermion and of the slave rotor, in 
the spirit of the non-crossing approximation. 



in which the average values are calculated with the ef- 
fective Hamiltonians above. In addition, the Lagrange 
multiplier h is determined from the constraint equation: 

We =E -5) (36) 

Let us emphasize that, in the decoupling leading to 
(|33I34|I . we have assumed that the average values 

(/il/jc) / an d (expi(^ — 9j))e on a given bond are both 
real. In fact, one could look for more general classes 
of solutions in which both (f- a fja - f- a fia) f ^ and 
(sin((?i — Oj))e 7^ 0. This would correspond to solutions 
with orbital currents around a plaquette, as proposed by 
several authors [25|. Spontaneous orbital currents are 
very naturally described using the slave rotor method, 
but will not be considered further in this paper, which 
aims at the general formalism. 



B. Hubbard model 

1. Slave rotor formulation 
In this section, we consider the Hubbard model: 

H = J2 H at(i)-J2 t H d lA<r (31) 
i ij,er 

which can be rewritten in terms of the rotor and auxiliary 
fermion variables as: 



U 



(32) 



Note that, in this context, — e = /1 is the chemical po- 
tential controlling the average density per site. Let us 
make a first approximation, which consists in decoupling 
the rotor and fermion variables on links (bes ides treating 
the constraint on average, as above), see [24| for a similar 
approach in the case of the t-J model. We then obtain 
two effective Hamiltonians: 

H f = - E 1 f fU* + ( e ° - h ) E ( 33 ) 



corresponding respectively to free fermionic spinons with 
an effective hopping and to a quantum XY- model 
for the phase variables with effective exchange constants 
. These effective parameters are determined by cou- 
pled self-consistent equations: 

tf = fy (cos^ - 9 S )) B , Jf = E *« (35) 



2. Simplest mean-field 

In the next section, we shall investigate some physi- 
cal consequences of equations (|33I34I35|I which approx- 
imate the Hubbard model by free spinons coupled self- 
consistently to an XY-model for the phase degrees of free- 
dom. We point out that the decoupling between fermion 
and rotor degrees of freedom can be viewed as a con- 
trolled approximation corresponding to a large-N limit 
of a multichannel model, as detailed in appendix [X] 

Here, in the same spirit as above, we consider a fur- 
ther simplification, which consists in treating the quan- 
tum XY model at the mean-field level. In this framework, 
the phase degrees of freedom is described by a mean-field 
Hamiltonian of independent sites: 



H, 



MF 



E 



— Li" 2 
2 



hLi + K cos 6i 



(37) 



with K = -2 J2j J if (cos 9j) e . Combining this with (jSHJ) 
and calculating the average values with the free-fermion 
Hamiltonian (|33|) . we finally obtain the following self- 
consistency equations for the variational parameters K 
and h: 

K = 2N(cos6) JdeD(e)en F (Ze + e - h) 



(L) = N deD(e) 



n,F (Ze + eo — h) 



(38) 
(39) 



Finally, the relation between the chemical potential — eo 
and average number of particle per site and color n is 
given by: 

n = Ytf'M =J de D ^ np ( Ze + 60 - V ( 4 °) 
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In these expressions, D(e) = J d-™ 6(e — e^) is the den- 
sity of states (d.o.s) of the band in the absence of in- 
teractions. The auxiliary fermion (quasiparticle) Green's 
function reads: 



Gf (k, icj n ) 1 = iuj n - e + h - Z e k 



(41) 



We recognize Z as the quasiparticle weight, which also 
determines the quasiparticle mass enhancement m* jm = 
\jZ. These two quantities are related because of the 
simple single-site approximation made here. 

At zero temperature, the number equation (|40|l implies 
that: 



h- e = Z no{n) 



(42) 



in which (Xq is the chemical potential of the non- 
interacting system such that J_ deD(e) = n. From 
(|41|l . it is seen that the Fermi surface is located at 
£k = (h — £o)/Z, and thus Q42|l implies that Luttinger 
theorem is satisfied. In fact, within this simple approx- 
imation in which the self-energy is independent of mo- 
mentum, the Fermi surface is unchanged by interactions 
altogether. The equations for K and h, at T — and for 
a given density, simplify into: 



K 



with e(n) = J^°^ n ' deD(e)e the average kinetic energy 
per electronic degree of freedom in the non-interacting 
model. 



N n 



1 



2Ne(n) (cos6)t 



(43) 
(44) 



3. Mott transition and orbital degeneracy 

We expect a Mott transition to occur at each com- 
mensurate filling n — Q/N (Q being an integer). This 
is associated with the vanishing of Z, and therefore 
the above equations can be analyzed analytically close 
to the transition (where Z is small) from a perturba- 
tive analysis in K, similar to the one performed in sec- 
tion ^^^] for the Anderson model in the Kondo regime. 
The ground-state wave function of H^ F in the insulat- 
ing phase (Z = K = 0) is e u ™ e with t n = N(n - 1/2). 
First-order perturbation theory in K yields: 



2^E 



(VA COS 01*/ 



2UK 



Ef 



U 2 - 4(U£ n + h) 2 



0(K 2 



(45) 



Since Z vanishes at the transition, but Ho(n) is finite, it 
follows from (|42|l that h = eo- For vanishing Z, the 
relation between eo and n is identical to that of the 



atomic limit, Eq. I|13(l established in the previous sec- 
tion: l n = Int [1/2 - e /U] with n = l/2+£ n /N. Finally, 
combining 145(1 and 144JI . we obtain: 

U c {e ) 2 - 4 [U c (e )£ n + e ] 2 + ANI (n)U c (e ) = (46) 

In this expression, t n and n should be viewed as depend- 
ing on the chemical potential e according to the relations 
just given. This expression determines the boundary 
Uc(co) between the metallic and Mott insulating phase in 
the (eo, U) plane. It is depicted for the case N = 4 (two 
orbitals with spin) in Fig. 01 The condition dU c /deo = 



U/D 




FIG. 3: Phase diagram for iV = 4 (two orbitals) at T = 
0, as a function of the chemical potential eo = —fj, and the 
interaction strength U/D. The three lobes correspond to the 
Mott insulator phases associated with half-filling (Q = 2) and 
quarter- filling (Q = 1, 3) respectively. 

determines the tip of each insulating lobe, i.e the criti- 
cal coupling U c {n) at which the insulating phase is en- 
tered as one increases U for a fixed commensurate density 
n. Differentiating (|46|) . it is seen that this happens for 
eo = —Ul n , i.e precisely at the center of each step of the 
Coulomb staircase. The critical coupling thus reads: 



U c (n) = W\e(n)\ 



(47) 



The phase diagram in the (n, U) plane is depicted in 
Fig. 0] for a flat d.o.s of half-width D, in which case 
U c = 4NDn(l — n). We see that the critical coupling 
is biggest at half-filling n = 1/2 (Q — N/2), which is ex- 
pected since orbital fluctuations are largest in this case. 
This conclusion may depend on the precise shape of the 
d.o.s however (and in particular may not hold for den- 
sities of states such that D(—e) ^ D(e)). The critical 
coupling increases linearly with orbital degeneracy N. In 
fact, an analysis of the DMFT equations for large orbital 
degeneracy was made in Ref. |ll|. and the exact behav- 
ior of the critical coupling at leading order in N found 
there is correctly reproduced by the simple mean-field de- 
tailed here. It is also instructive to compare the present 
results with that of the multi-orbital Gutzwiller approx- 
imation [26|, which reads: J7 C GA = i(N + 2)|e(n)|. Our 
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FIG. 4: Phase diagram in the (n,U) plane. The Mott insula- 
tor lobes collapse to lines at commensurate fillings, when U 
is larger than U c {n) (shown as dots). 



expression has the same behavior at large N, but yields in 
general a smaller critical coupling: U c — U^ A N/ (iV + 2). 
For small orbital degeneracies, we believe (on the basis 
of, e.g, DMFT results) the Gutzwiller expression of U c 
to be quantitatively more accurate. 

The slave-rotor mean field equations are easily solved 
numerically by determining iteratively the parameters h 
and K. At each iteration, the spectrum of the single- 
rotor Schrodinger equations is computed (using e.g a de- 
composition on the atomic basis states e lle ). In Fig. EJ 
the ground-state wave function ^o(8) is displayed for sev- 
eral values of U at half-filling. The curves nicely illustrate 




FIG. 5: Rotor ground state wave function \I/o(#) with values 
of the local interaction ranging from U/U c = 0.01 (peaked 
curve) to U/U c — 1 at the Mott transition (flat curve). 

how one goes from the insulator (in which case there are 
little charge fluctuations, and maximal phase fluctuations 
so that the wave-function is delocalized over all 9 values) 
to the metal (in which case charge fluctuations become 
large at small U, and the wave-function is peaked such 



as to limit phase fluctuations) . The corresponding quasi- 
particle weight is displayed in Fig. as a function of 
U/Uc- The simple slave-rotor mean-field is compared to 
the DMFT result and to the Gutzwiller approximation 
(GA). It is seen that, close to the transition, the slave- 
rotor mean field reproduces more accurately the DMFT 
answer than the GA. It is not very accurate at weak- 
coupling however (even though Z correctly goes to Z = 1 
at U — 0, it has an incorrect small- U expansion). In fact, 
it is a quite general feature of this slave-rotor mean field 
that the method is more accurate in strongly correlated 
regimes. 
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FIG. 6: Quasi-particle weight Z as a function of U/U c at 
T = 0; DMFT calculation (thin line), rotor mean-field theory 
(thick line) and Gutzwiller approximation (broken line). 

In Fig. [3 we plot the number of particles as a func- 
tion of the chemical potential for N — 4. The value 
of U has been chosen to be bigger than the critical cou- 
plings yielding an insulating state, for any commensurate 
filling. The curve illustrates the plateaus found at each 
commensurate filling, the central one (half-filling) being 
narrower (compare to Fig. The effective mass en- 
hancement (= 1/Z) is also plotted in Fig.|Hlas a function 
of chemical potential for a smaller value of U, such that 
a metallic phase is found at any filling. The curves illus- 
trates how a largest effective mass enhancement is found 
at low and high fillings n = 1/4,3/4, and a compara- 
tively smaller close to half- filling n — 1/2 (again, this 
conclusion depends on the shape of the d.o.s). 

The description of the Mott transition obtained within 
this simplest mean-field has many common features with 
the Brinkman-Rice (BR) [2]} one. Indeed, the effective 
mass diverges at the transition and the quasi-particle 
residue vanishes [Z ~ 1 — 17 /U c ) as in BR. There is 
one significant difference however, which is that in the 
present description, the optical gap A of the insula- 
tor does not coincide with the chemical potential jump 
A/i = — Ae = ji{n + + ) — /i(n — + ) for infinitesimal 
doping away from a commensurate filling. Indeed, within 
this simple mean-field, the spectral function of the insu- 



9 




FIG. 7: Total occupancy Q = 4n as a function of eo for a 
value U = 4.5 larger than all critical interactions U c (n), in 
the two orbital case (N — 4). The Mott insulators are seen 
here as charge plateaus. 
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FIG. 8: Effective mass m*/m for U = 2.5 below all Mott 
transitions, in the two orbital case (N — 4). 

lator is identical to that of the atomic limit (not sur- 
prisingly, the simple mean-field with only two variational 
parameters describes the charge fluctuations in the in- 
sulator in an oversimplified manner). As a result, the 
optical gap simply reads 

A = U (48) 

in our approach, and is therefore not critical at the Mott 
transition. In contrast, the chemical potential jump van- 
ishes continuously at U c . Indeed, solving l|46|l for eo 
yields: 



Am = ^\/1- 77 ( 49 ) 



These features are very similar to those obtained within 
dynamical mean-field theory [2|. This is not surprising, 
since the single-site mean field approximation to the XY- 
model indeed becomes exact in the limit of infinite coor- 
dination of the lattice. Note however that this is not the 
case of the approximation (|33|I - H^5(I which consists in de- 
coupling the rotor and fermion variables (see section ll*V|l . 
Within DMFT, the quasiparticle weight vanishes at a 
Brinkman-Rice like critical point U C 2 while the optical 
gap of the insulator vanishes at a Hubbard-like critical 
point U c \. As a result, the strongly correlated metal 
close to the transition displays a clear separation of en- 
ergy scales: the quasiparticle coherence scale e* F ~ ZD 
being much smaller than the ("preformed") gap of the 
insulator A. The simple mean-field of this section is in a 
sense a somewhat extreme simplification of this picture, 
in which U C 2 = U c and U c \ is sent to U c \ = (this is 
consistent with the known fact [HI that U c2 cx N while 
U c i oc VA, and that the simple mean-field becomes more 
accurate for large- A). 

IV. INCLUDING SPATIAL CORRELATIONS 
AND PHASE FLUCTUATIONS 

In this section, we go beyond the single-site mean- 
field approximation, and investigate the physical conse- 
quences of the approximate description of the Hubbard 
model introduced in Sec. IIII B II This description, sum- 
marized by Ecis. l|33l35|) . consists in a free fermion model 
Hf coupled self-consistently to a quantum XY-model Hg 
for the phase degrees of freedom. 

A. General considerations 

Let us first emphasize some general aspects of this de- 
scription, before turning to explicit calculations. The 
Hamiltonian for the phase degrees of freedom has two 
possible phases: a disordered phase without long-range 
phase order, and a long-range ordered phase. At zero- 
temperature, one expects a quantum phase transition 
from the ordered phase to the disordered phase as the 
ratio /7/j7jj is increased. Since the Green's function of 
the physical electrons read, within this approximation: 

GUr-r') = G{ j (r-r')(e^- s ^) e (50) 

it is seen that the quasiparticle weight Z, associated with 
the limit of large-distance and large time separation (low 
frequency), is given by: 

Z = (cosO^l (51) 

Thus, the phase with long-range order for the rotors cor- 
responds to the metal (Z ^ 0), while the disordered 
phase correspond to the Mott insulator [Z = 0). Ob- 
viously, the description of the Mott metal-insulator tran- 
sition that follows is closely analogous to that of the 
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superfluid-Mott insulator transition in the bosonic Hub- 
bard model 28, 29] . Two remarks about this description 
of the metal and of the insulator are in order. First, it is 
of course unphysical to think of a metal as having long- 
range phase coherence. Naturally, this is only true of 
the saddle-point approximation in which the rotors and 
spinon degrees of freedom are entirely decoupled. Fluc- 
tuations will induce interactions between these degrees 
of freedom, restore inelastic scattering and thus destroy 
phase coherence. The absence of inelastic scattering at 
the saddle-point level is a well-known feature of slave- 
boson theories. Note futhermore that despite the order- 
ing of the rotors, the metallic phase becomes a supercon- 
ductor only when (//t/Ji ) is also non-zero (i.e when there 
is spinon pairing). Second, the insulator envisioned here 
is a non-magnetic insulator without any spin or trans- 
lational symmetry breaking, i.e a spin-liquid. Even in 
the disordered phase, (cos(#, — 0j)) on a given bond (e.g 
nearest-neighbour) is non-zero (it corresponds to the en- 
ergy density of the XY model). Therefore t\^ ^ in the 
insulating phase, so that the spinons have a Fermi sur- 
face (with Luttinger volume). This also implies that t\^ 
remains finite through the Mott transition and therefore 
that the effective mass does not diverge, despite the fact 
that Z — > 0. These last remarks apply to any finite di- 
mension, but of course not to d = oo. In this limit, the 
single-site mean field of the previous section applies and 
(exp(9i — 9j)) = (cos#) 2 = Z. Finally, we emphasize that 
the non-magnetic nature of the insulator is of course as- 
sociated with the fact that the rotor degrees of freedom 
are associated with the charge and are not appropriate to 
properly describe spin ordering. Therefore, they are bet- 
ter suited to lattices with strong frustration (or models 
with large orbital degeneracy) in which a spin-liquid insu- 
lator is a realistic possibility. Finally, because long-range 
order for the rotors corresponds to breaking a continu- 
ous 0(2) symmetry, a Goldstone mode will be present 
in the ordered (metallic) phase. This mode is present in 
any finite dimension, but disappears in the d = oo limit. 
It corresponds to the zero-sound mode of the metal. As 
we shall see, these long-wavelength modes play an im- 
portant role: they change the low-energy description of 
the transition as compared to the d — oo (DMFT) limit. 
As a result, the separation of energy scales does not ap- 
ply in a strict sense (the "preformed" gap found within 
DMFT is filled up with spectral weight coming from these 
low-energy modes). As we shall see however, this spec- 
tral weight remains small in high dimensions (including 
d = 3), so that an approximate separation of scales still 
applies. 



B. Sigma-model representation: saddle-point 
equations in the spherical limit 

In order to perform explicit calculations with the quan- 
tum rotor Hamiltonian (|34fl . we shall use an approxima- 
tion that has proven successful in the context of quantum 



impurity models with slave rotors [7], |8| • It consists in re- 
placing the quantum rotor exp(i#i) by a complex bosonic 
field Xi(r) and to treat the constraint \Xi\ 2 = 1 on aver- 
age. Alternatively, this can be viewed as extending the 
0(2) symmetry to O(M) and taking the large-M (spher- 
ical) limit. This is a well known approximation to non- 
linear sigma models 29], which preserves many qualita- 
tive features of the quantum phase transition. For details 
of the formalism in the slave rotor context, see Ref. Q. 
In the following, we focus on the half-filled case (since 
we are mainly interested in the Mott transition), with 
a particle- hole symmetric d.o.s D(e), so that we can set 
en = h = 

The spinon and rotor (now X-field) Green's function 
read |39| : 

G/(k, iujn)" 1 = iuj n - Q/e k (52) 
Gx(k,z^„)- 1 = ^ + \ + Q x e k (53) 

In these expressions, w„ and v n are respectively, 
fermionic and bosonic Matsubara frequencies, A is a La- 
grange multiplier associated with the constraint (|AT| 2 ) = 
1, while Qf and Qx are the self-consistent parame- 
ters entering the effective spinon hopping and XY cou- 
pling constants: Qf — (cos(#; — 9j)) = (XiX*), Qx — 

fiafiv) • The self-consistent equations which deter- 
mine A, Qf and Qx read: 

1= / deD(e)-Y^- — (54) 

J-D vl/U + X + Qxt 

DQf = - [ deD(e)e^Yl 2 /TT A . n (55) 
J-d P „ K/ u + A + Qxe 

DQ X = -2 f deD(e) en F {Q f e) (56) 

J—D 

These expressions have been written here for a simple 
tight-binding band with nearest-neighbor hopping ty = t 
on a d-dimensional cubic lattice (e^ = —It Y] a —i cosk a ). 
As above, D(e) denotes the band d.o.s, and D = 2dt 
is the half-bandwidth. For simplicity, we have set the 
orbital degeneracy to N = 2 in these equations. 



C. The Mott transition: Mott-Hubbard meets 
Brinkman-Rice 

In this section, we investigate the solution of these 
equations at zero-temperature. This leads to a descrip- 
tion of the finite-dimensional Mott transition that we an- 
alyze in detail. 
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1. The insulating phase 

Lets us note first that Eq. (|56*)l readily determines Qx 
at T = 0: 

DQx | T=0 = -2 / £>(e) e de = 2|e | (57) 

•/-£> 

From the form i|53|) of the X-field Green's function, one 
sees that the bosonic spectrum has a gap as long as A — 
Qx D = A — 2|e| > 0. In this case, there is no long- 
range order for the phase degree of freedom, and this 
corresponds to the insulating phase. The insulating gap 
reads: 



A 3 = 2y/U(\-Q x D) 



(58) 



and we can rewrite Eqs. I|54I55|I as selfconsistent equa- 
tions for the gap A g and the renormalization of the 
spinon hopping Qf. This reads, at T = 0: 



1 = f (ieD(e) 

J-D 

Qf = - f°de D(c) 

J-D 



U 



A 2 +8U\e\(l + e/D) 



c 

D 



U 



^l + W\e |(1 + e/D) 



(59) 



(60) 



These equations are valid in the insulating phase, when 
A 9 > 0. The gap vanishes at a critical coupling U c ob- 
tained by setting A g = in 



Uc 



D 

de _ 

D y/T 



D(e) 



1 -2 



e/D 



(61) 



In this expression, U£° — 8|e| is the critical coupling 
corresponding to the d = oo limit, in agreement with ex- 
pression (|47|l of the previous section (with N = 2). Note 
that, in the d — > oo limit, one must scale the hopping as 
t = t* /Vd, so that D = 2dt oc \fd — > oo and the r.h.s of 
(|61|l goes to unity. The integral in l|61(l is smaller than 
unity in general, so that U c decreases as dimensionality 
is reduced. We also note that in one dimension, this in- 
tegral has a logarithmic singularity at band edge, since 
7rD(e) = D/VD 2 - e 2 , so that Eq. ^ yields C/ C d=1 = 0, 
which is indeed the exact result for a half-filled Hubbard 
model with N = 2 [3(j (see however jiof'). 

Substracting Eq. I|59|l from the same equation with 
A g = (which defines C/ c ), one obtains: 



D 



u^_ 
u 



(62) 



l-D 



de D(e) 



V 1 + e / D ^A 2 g /(UU^) + l + e/D 



The expansion of this expression for small A g depends 
dimensionality. For d > 3, the integral J deD(e)(l + 



on 



f-/D) 3 / 2 is convergent at band edge e = —D, recalling 
that D(e) ~ (D + e) d / 2 " 1 near the bottom of the band. 
In contrast, the small-Ag expansion is singular for d < 3. 
This analysis finally leads to the following behavior of 
the gap close to the critical point: 



A g /Uc oc 



u 



1 for d > 3 



i/(d-i) 



for d < 3 (63) 



Hence we find that the exponent changes from its mean- 
field value 1/2 for d > 3 (as found e.g in the single-site 
mean-field of the previous section and the Gutzwiller ap- 
proximation) to a non- mean field exponent for 1 < d < 3. 
Therefore, d — 3 corresponds to the upper critical dimen- 
sion in this description of the Mott transition (logarith- 
mic corrections are found in that case). Below d = 3, the 
exponent l/(d— 1) corresponds to that of the large-M 
limit of the quantum O(M) model in d-dimensions, i.e 
to that of the d + 1-dimensional classical model. Had 
we kept 0(2) quantum rotors, we would have found 
Ag ~ (U /U c — l) zl/ with z = 1 and v the correlation- 
length exponent of the d + 1-dimensional classical XY 
model, as in the case of the superfluid-insulator transi- 
tion of the Bose Hubbard model 281. 



2. The metallic phase 

For U < U c , the gap closes and one enters the metallic 
phase. In this regime, the constraint equation (|X| 2 ) = 
1 can only be satisfied by a Bose condensation of the 
X-field. As in studies of quantum magnetism based on 
Schwinger bosons |3l| , Bose condensation in the spherical 
limit corresponds to the phase with long-range order for 
the rotors. In this phase, the constraint equation l|54|) 
has to be rewritten by isolating the k = mode in the 
Brillouin zone. The Lagrange multiplier A sticks to the 
value A = QxD — 2|e| in this phase. The full X-field 
Green's function thus reads at T = 0: 



G x (k,iv) - ZS(v)6{k) 



uyU + 2\e\(l + e k /D) 



(64) 



The condensation amplitude Z — (X) 2 is determined 
from the constraint (\X\ 2 ) = £k G *( k > T = 0) = 1: 



U 



-D yoe(l + e/D) 

which simply reads, using H61(l : 

IT 



Z = l -^U C 



(65) 



(66) 



This expression vanishes linearly, Z ~ (U c — U)/2U C , at 
the critical point for all dimensions d > 1. The fact that 
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there is no change of critical behavior for Z at d = 3, 
in contrast to the gap, is due to the use of the spherical 
approximation. Had we kept 0(2) rotors, we would find 
Z ~ (1 - U/U c ) v ( d - 1+r >), with v and n the critical expo- 
nents of the d+ 1-dimensional classical XY model. In the 
spherical approximation rj = and v = 1 / (d — 1 ) so that 
Z ~ (1 - C//C/ c ) also below d = 3. 

Hence, we have found that the quasiparticle weight 
and insulating gap vanish at a unique critical coupling. 
As shown below, the gap A g given by is the gap 
in the single particle spectral density of the insulator. It 
also coincides with the chemical potential jump Afi when 
the present approach is extended away from half-filling. 
Hence, in this finite dimensional description of the Mott 
transition, we find a unique critical point correspond- 
ing both to Brinkman-Rice physics |2jj (vanishing of Z) 
and to Mott- Hubbard physics |32| (gap opening). This 
is in strong contrast to the d = oo single-site mean field 
investigated in the previous section, and to the DMFT 
picture Q. Below, we show that this is due to long wave- 
length collective modes filling in the preformed gap, and 
investigate in detail how the previous picture is recovered 
in the (singular) d = oo limit. 

The equation (|55(l for the renormalization Qf of the 
effective hopping must be rewritten in the metallic phase 
to take into account the Bose condensed fraction. At 
T = 0, it reads: 

Qf = — = Z~ \U~ Tde D{e)^= (67) 

This expression makes very clear that the effective mass 
remains finite at the critical point, while Z — > (note 
that the integral in the r.h.s of l|67l) is negative so that 
Qf > Z). In the d — > oo limit, one recovers m/m* = Z, 
since D oc yd — > oo. This calculation can be extended 
to the weakly doped Mott insulator at large U and hole 
density S, with the result: 



to t/U + 5 J/t + S v ; 

Hence, the present theory correctly captures the mag- 
netic exchange energy J oc t 2 /U, which quenches out the 
spin entropy (due to spinon degrees of freedom) in the 
insulator and hence prevents the effective mass from di- 
verging at the Mott transition. This is expected from the 
fact that the spinons form a dispersive band in the insu- 
lating state and thus have an entropy depending linearly 
on T at low temperature. These findings are entirely 
similar to the slave bosons mean-field theories of the 
t— J model ;33]. Fig. [5] and Fig.llUI illustrate graphically 
the physical quantities characterizing the Mott transition 
which we discussed previously. 




U/U c 



FIG. 9: Plot of the quasiparticle weight Z, the effective mass 
renormalization Qf = m/m* and the Mott gap A 9 as a 
function of U /U c across the Mott transition in the three- 
dimensional case. 
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FIG. 10: Effective mass m* /m = 1/Q provided by the mean- 
field equations I|54I56|I for d — 3 (bold line) and d — oo (thin 
line). For comparison, a DMFT-IPT calculation (dashed line) 
is also presented. 



3. Spectral functions and collective modes: what are 
Hubbard bands made of. 

The Green's function of the physical electron in the 
approximation of decoupled spinons and rotors is given 
by {HJ as: Gf^r) = G^^G^ir). Using JH1 for the 
X-field Green's function, this leads to: 

G d (k, iw) = , f- + GT(K iw) (69) 

IUJ - Qf £k 

This expression is valid in the metallic phase. In Fig. II II 
we display the k-integrated (local) spectral function 
^ k Pd(k, cj), as the Mott transition is approached, in 
the three-dimensional case. The first term in l|69|l cor- 
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1.5- 




FIG. 11: Zero temperature local density of states across the 
Mott transition for a three-dimensional cubic lattice, with 
[7 = 0, £4/2, U c , 3C/c/2 (for D = 1). 



responds to the coherent quasi-particle. When summed 
over k, the quasiparticle contribution to the local spec- 
tral function yields a peak Z/QfD(u)/Qf). The spectral 
weight of this peak is Z, its width is of order QfD and 
its height is ZD(0)/Qf. Hence, its height goes to zero 
at the transition, while its width is reduced but remains 
finite (Fig. Ill|) . In the d = oo limit (where Z = Qt) 
the zero-frequency density of states is pinned at its non- 
interacting value, as known from the Brinkman-Rice pic- 
ture and the dynamical mean-field theory: only in infi- 
nite dimension does the quasiparticle peak disappear by 
narrowing down instead of collapsing. 

The incoherent contribution G™ c comes from the con- 
volution of the free spinon Green's function with the 
non-condensed contribution to the rotor Green's func- 
tion, G™ c = l/y/U + 2|e |(1 + eu/D)], which is the 
second term in l)64[) . The latter corresponds to bosonic 
collective modes dispersing according to: 

The incoherent contribution corresponds to the Hubbard 
band, which are well developed in the correlated metal, as 
also predicted by DMFT. Note that the dispersing branch 
of bosonic excitations is centered around ±WUU°° /2. 
Hence a measure of the typical energy scale associated 
with the distance between the two Hubbard bands is: 

A = ^fUXjf (71) 

However, a key point is that this branch of collective 
modes extends to arbitrary low frequency where it be- 
comes the Goldstone mode of the broken symmetry. In 
the small momentum limit, the dispersion relation 171)11 
reads: 

luu°° 

Wx(k) ~ V"8d" |k| (k " 0) (?2) 



The corresponding density of states behaves as px{u) ~ 
u! d ^ 1 . These long wavelength excitations are responsible 
for tails of the Hubbard bands, extending down to low- fre- 
quency. This low-energy spectral weight due to collective 
modes is the origin of the continuous closure of the Mott 
gap at the Brinkman-Rice transition. In other words, the 
Hubbard bands are made of two kinds of contributions. 
The main part of their spectral weight is associated with 
bosonic modes whose momentum is not small, so that 
ujx is finite (of the order of A). In addition, the small 
weight in the (~ tails at low-energy is associated 

with the k ~ collective modes. In the Mott insulator, 
all the bosonic modes are gapped, but the bottom of the 
bosonic density of states has (w — Ag) d_1 tails which con- 
tribute to the tails of the Hubbard bands, filling in the 
energy range between A g and A. 

It is natural to interpret the bosonic collective modes as 
the zero-sound mode of the metal. Indeed, these modes 
have been discussed previously by Castellani et al. |33|. 
in their study of fluctuations around the saddle-point 
of conventional slave boson approaches. These authors 
pointed out that the Mott transition is associated with 
the softening of this collective mode, as also found here. 
In the present approach, the collective modes appear on 
the same footing than the quasiparticles. 



4- The d — oo limit, separation of energy scales and 
"preformed" Mott gap 

It is instructive to understand more precisely what 
happens as the dimensionality is increased. As clear from 
the previous discussion, the d = oo limit is singular in 
at least this respect that the long-wavelength collective 
modes are absent. Indeed the sound velocity in (|72|l van- 
ishes in this limit. In fact the bosonic modes no longer 
have a dispersion: k-dependence disappears from the dis- 
persion relation l|7U|) since D must be scaled as D oc yd. 
The bosonic spectral functions thus has two poles on top 
of the condensed fraction, which leads, after perform- 
ing the convolutions, to the following simple form of the 
physical k-summed local spectral density (using also that 
Qf = Z in this limit): 

(73)' 

In this expression, A is given by l|71|> and corresponds 
to the typical separation between the Hubbard bands. It 
is sometimes referred in the framework of DMFT as the 
"preformed gap" in the metallic state (i.e the Hubbard 
bands are well separated from the central quasiparticle 
peak). A does not vanish at the Brinkman-Rice point 
and beyond this coupling the insulator sets in with a 
finite gap A. Note that, in the present approximation 
where spinons and rotors have been decoupled, one sim- 
ply has A — U in the insulator (as found also in (|48l) ') 
and that, accordingly, the Hubbard bands in Ij73[) have 
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vanishing width close to the transition. This pathological 
result can be improved by including dynamical fluctua- 
tions of the auxiliary particles, as shown in Q . Despite 
these oversimplifications, the present approach does re- 
tain the main qualitative feature of DMFT, namely the 
separation of energy scales at the Mott transition. 

In Fig. 1121 we show how the large-d limit is approached 
by plotting the local spectral density right at the critical 
coupling U = U c , for increasing dimensionality. This plot 
clearly reveals the two-components building up the Hub- 
bard bands, with the main part of the spectral weight 
centered around the "preformed" gap A and tails ex- 
tending down to low-frequency (down to the true gap 
A g in the insulator), associated with the long- wavelength 
collective mode. The inset demonstrates that as dimen- 
sionality increases, the spectral weight in the tails be- 
comes smaller, so that an approximate separation of en- 
ergy scales holds (and in fact already holds to a good 
approximation in d = 3, while it is no longer meaningful 
in two dimensions). 




FIG. 12: Local spectral density at the critical point U = U c , 
for increasing dimensionality (top to bottom curves at small 
uj: d = 2,3,4,5). The spectral weight associated with the 
low-frequency tails of the Hubbard bands is seen to decrease 
as dimensionality increases. Correspondingly, the separation 
of energy scales and the pre-formed gap become more and 
more apparent. The inset shows the integrated density of 
states J Q depd(e) demonstrating this approximate separation 
of scales for d > 3. (Note that the progressive narrowing 
of the main lobe of the Hubbard band as d increases is an 
artefact of the approximation in which spinons and rotors are 
decoupled. 



V. CONCLUSION AND PERSPECTIVES 

In this paper, we have used the slave-rotor represen- 
tation in order to construct approximation schemes for 
strongly correlated fermion models. A theory of the 
Hubbard model involving free fcrmionic spinons self- 
consistently coupled to a quantum XY model has been 
developed. The quantum phase transition of the latter 
corresponds to the Mott transition between a strongly 



correlated metal and a Mott insulating spin-liquid with 
a spinon Fermi surface. Both the gap in the spectral 
function and the quasiparticle weight vanish at the criti- 
cal point, while the effective mass remains finite. In this 
picture, long-wavelength collective modes of the phase 
variable play an important role. They are responsible 
for low-energy tails of the Hubbard bands (in addition to 
the main component of these bands, which are atomic- 
like short-distance excitations). This has potential impli- 
cations for spectroscopic and tunneling experiments. In 
infinite dimensions, these collective modes arc suppressed 
and this limit appears singular in this respect. Only in 
this limit is a strict separation of energy scales recovered, 
as in the DMFT picture. 

This raises some open questions associated with the 
physics of these collective modes, which are physically 
associated with zero-sound. At the saddle-point level, 
where spinons and rotors do not interact, these modes are 
undamped. The metallic state is described as a perfect 
metal with no inelastic scattering. Interactions between 
rotors and fermionic spinons will induce a Landau damp- 
ing of these modes, and the metallic state will loose phase 
coherence. One possibility is that this damping is large, 
which would presumably weaken its effect and might re- 
store some separation of energy scales as in the DMFT 
picture. This calls for further work on the nature of col- 
lective modes in a metal close to the Mott transition. In 
particular, the restauration of gauge symmetry, broken 
at the saddle-point level, will have a strong impact on 
the non-gauge invariant pr opag ators, as already known 
for single impurity models |34|. However, it is yet not 
entirely clear whether this effect will affect strongly the 
physical sector at low energy. One additional issue is that 
in a real material the screening deteriorates as one gets 
closer to the Mott insulating state [3!|. The "acoustic 
plasmon" mode of the Hubbard model with short-range 
interactions will be pushed to higher energy and this may 
also weaken its relevance for low-energy physics. Another 
issue which will arise when taking into account the inter- 
actions between spinons and rotors is the description of 
the insulator as a non-magnetic spin-liquid. Stabilizing 
such a spin-liquid state beyond saddle-point level is pre- 
sumably possible only on a very frustrated lattice. Even 
in this case however, it has been suggested recently [3(| 
that a superconducting phase can intercalate between the 
metal and the insulator, due to the proliferation of short 
range spin singlets, therefore superseding the zero tem- 
perature Mott transition. Nevertheless, some of our re- 
sults, such as the finiteness of the effective mass, should 
remain valid above the low temperature ordered regions. 

Finally, we point out that the slave rotor representa- 
tion explored here is a useful technical tool that can be 
applied to strongly correlated systems in a variety of con- 
texts. Because a single collective variable is introduced, 
(which has a direct physical interpretation in connection 
with the local charge), using this representation is gen- 
erally simpler than other finite-U slave-boson schemes 
provided one deals with a symmetric interaction. Ap- 
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plications to mesoscopic devices and quantum impurity 
models have been presented elsewhere @, ■ Other po- 
tential applications are the effect of long-range or time- 
dependent interactions [22^. or the interplay of disorder 
and interactions. Interacting boson models can also be 
expressed with slave rotors (for a recent application of 
variational approximations to the XY model, in the con- 
text of bosonic models see 37]). Although mean-field 
approximations for interacting bosons can be formulated 
in a simple manner due to the commuting nature of the 
physical degrees of freedom [2j|, the slave- rotor repre- 
sentation might prove useful in the context of interacting 
cold atoms in order to deal, for example, with boson- 
fermion mixtures. 



APPENDIX A: LARGE-N LIMITS AND 
MEAN-FIELD APPROXIMATIONS 

We discuss here how the different mean field approx- 
imations presented in this paper can be formulated in 
terms of large N limits of generalized Hubbard models. 
Let us introduce a "multichannel" version of the Hub- 
bard model, based on spin-carrying fermions /J (a — 
1, • • • , N) and channel-carrying phases 8 a (a = 1 . . . K): 

icr ia ijcra 

(Al) 

Two Hubbard- Stratonovich fields conjugate to f^ a fja 
and J2 a e l ^ 9ia ~ e ^ can be introduced in order to decou- 
ple the last term. When both N and K are large, with 
a fixed ratio K/N, a saddle-point applies which leads to 
the decoupled effective Hamiltonians 1)3313411 . This cor- 
responds to a factorization the hopping term on bonds, 
as shown by the effective parameters (|35ll . A similar re- 



mark applies in the usual context of slave-bosons for the 
t — J model: the mean-field approximation investigated 
e.g in 24] correspond to a multichannel limit of: 

Ujfiafjgbiabja (A2) 

ijaa 

Because the quantum XY model on the lattice is not 
easily investigated analytically, we have performed in sec- 
tion ^^B] a 0(2M) generalization of the phase part that 
leads to further simplications, while allowing to deal with 
the model in finite dimensions. This can also be seen as 
a direct large N,M limit of a SU(A)xO(2M) Hubbard 
model (see |7j for a related approximation concerning the 
Anderson model): 

H = E f 0/il/i<r+^: E (^i" 3 ) _ E M /il/jff X *a X ja 
icr ia/3 ij^a 

.(A3) 

introducing a complex field with M colors X- la subjected 
to the spherical constraint: J2 a |Ai Q | 2 = M. In the pre- 
vious expression, L" 13 denotes the 0(2M) angular mo- 
mentum tensor associated with the X- la vector |29j . 

Finally, we note that the simplest single-site mean field 
of section IIII B 21 can also be seen as a large-N limit of a 
generalized Hubbard model which reads (note the differ- 
ent indices position and the scaling of the hopping term): 



icr ia ijaaa' 

. (A4) 

This gives a on site factorization of the phase variables. 
Alternatively, this can be seen as a large connectivity of 
the bond mean field approximation (|33I34() . 
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